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Abstract 

We compute the leading and next-to-leading corrections to the finite temper- 
ature scalar potential for a 3+1 dimensional (/>^ theory using a systematic 
expansion. Our approach automatically avoids problems associated with infrared 
divergences in ordinary perturbation theory in h. The leading order result does not 
admit a first order phase transition. The subleading result shows that the exact 
theory can admit at best only a very weak first order phase transition. For = 4 
and weak scalar coupling we find that Ti, the temperature at which tunneling from 
the origin may begin in the case of a first order transition, must be less than about 
0.5 percent larger than T2, the temperature at which the origin changes from being 
a local minimum to being a local maximum. We compare our results to the effective 
potential found from a sum of daisy graphs. 



1. Introduction. 



There has been much recent interest in the nature of the electroweak phase transition 
motivated by the possibihty of baryogenesis within the standard model itself. As noted 
by Kirzhnits amd Linde [1], Weinberg [2] and Dolan and Jackiw [3] over 18 years ago, 
a spontaneously broken field theory may have its symmetry restored at high enough 
temperature. For example, in a spontaneously broken 0^ theory, temperature corrections 
give a positive mass-squared contribution at the origin (0 = 0), and at high enough 
temperature this correction results in a global symmetry-unbroken minimum at the origin. 

In such a model there are two important temperatures as we lower the temperature 
from a very high value. The first, Ti, is the temperature at which a second possible 
minimum appears degenerate in energy with the minimum at the origin. The second, T2, 
is the temperature at which the effective mass at the origin vanishes, i.e. when the origin 
changes from being a local minimum to a local maximum. 

The phase transition from the symmetric to nonsymmetric phase as we cool a system 
described by such a model may therefore proceed in two ways, by tunneling when the 
temperature is between Ti and T2, or by a rollover when the temperature drops below 
T2. At T = Ti, isolated bubbles of the symmetry broken phase will start to be created 
by tunneling. If the phase transition completes by bubble nucleation we will call the 
transition first order, otherwise if it completes mostly by a rollover after T < T2 then we 
will call the transition second order. If the finite temperature scalar potential possesses 
only the local minimum at 0^ = until T = T2 then there is no Ti and the phase transition 
is necessarily second order. Even if there is a Ti we will say only that the system admits 
a first order transition since if Ti is close enough to T2 we expect the phase transition to 
complete by a rollover. 

In the context of the standard model, it has been suggested that sufficient baryoge- 
nesis may occur [4,5,6] if the phase transition is first order. Unfortunately, an accurate 
determination of the the nature of the phase transition has proven difficult in this case. 
One needs a reliable determination of the temperature dependent effective potential of 
the scalars near the origin for temperatures between Ti and T2. One-loop calculations 
suggest that baryogenesis can only occur in the minimal standard model for Higgs mass 
Mff < 45 GeV [7,8], a possibihty that appears ruled out by experiment [9]. However, as 
noted by both Weinberg [2] and Dolan and Jackiw [3], naive finite temperature pertur- 
bation theory for even the four dimensional 0^ model suffers from infrared divergences. 
These divergences invalidate the one-loop approximation. In particular, at temperatures 
close to T2 perturbation theory in h begins to break down. The reason for this can be eas- 
ily understood: even the one-loop correction drastically modifies the tree level potential; 
it is clearly not a small perturbation. Due to such problems, the one-loop result in pure 

theory gives a complex potential for small values of 0^ and cannot be used to study 
the nature of the phase transition. 

In the standard model, one may argue that if the phase transition completes "well 
before" T2 then the one-loop calculation can be used. However, in the case of the standard 
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model the phase transition is expected to be weakly first order, if it is first order at all. In 
addition, it is not a priori clear that resumming some infinite class of diagrams does not 
rule out a first order transition. The task then is to extract the leading corrections near 
T — T2 io all orders in perturbation theory. For example, Dolan and Jackiw summed a 
class of diagrams, the "super-daisies" to circumvent the infrared divergence problem in 
pure 0^ theory (with N scalars) to get a reliable estimate of T2. Such diagrams correspond 
to iterating the daisy graphs 




(1) 

so that each petal has an arbitrary number of petals and so on. Such a summation 
corresponds to the leading term in a 1/N expansion and greatly modifies the 1-loop 
estimate. 

In this article we will use the systematic expansion method to calculate the leading 
and next-to-leading corrections to the finite temperature effective scalar potential near 
the origin in a weakly coupled 0^ theory in four dimensions. In the case of the standard 
model, one also needs to incorporate corrections to the scalar potential from gauge loops. 
We do not address this problem here, but will later comment on the apphcability of 
our methodology to a gauged 0^ theory. In section two we review the 1/N expansion 
for a 0^ theory and present our computations. In section three we discuss our results. 
Qualitatively, we find that the leading order in N correction does not admit a first order 
transition. The next-to-leading may admit a first order transition, but an extremely weak 
one. We stress that this approach automatically avoids the infrared divergence problem 
of ordinary perturbation theory. In particular, for the range of (fP' and T we are interested 
in the potential does not become complex. 

In the context of the standard model there have been several recent papers on improv- 
ing the 1-loop estimate [10,11,12]. These attempts mainly involve summing the daisy 
diagrams. The simplest procedure [13] has been to use the 1-loop generated temperature 
dependent mass to calculate 1-loop corrections. In pure (f)^ theory, if properly done, this 
can be interpreted as summing the daisy diagrams. This procedure has its roots in the 
work of Weinberg [2] . It is useful to briefly discuss this procedure, and at the same time 
introduce some necessary formulas, before presenting our results. 

To understand the meaning of such a procedure, it is necessary to recall an important 
formula in the background field method. To be precise we consider a simple 0{N) scalar 
field theory with action {i = 1 . . . N) 

= \j s^^d^d^cp' - 1 /('^^ - ^')'- (2) 
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One calculates [14] 



m 



-zhln I m exp ('-{S[cl> + 0] - ^00^ 



h 



6(1)' 



'^'-^(m)- (3) 



The first term on the RHS of the last line is the classical action. The inverse scalar 
propogator is 

A^ = [M' + MJ(0)]. (4) 
Mf, is the second derivative of the classical potential w.r.t. the fields 0*, M^, — didjV — 



Equation (3) is an effective action in that it incorporates the effects of all one-particle- 
irreducible (IPI) diagrams and can be given to any order in perturbation theory. It may 
be evaluated in perturbation theory as 



(5) 



The second term is the familiar one-loop contribution and the rest contains all higher 
loop corrections from IPI graphs. The one-loop corrected potential is 
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V{<P) = V{(j)) - -i I Trln - 



(6) 



where the trace is over internal indices and for zero temperature the integration measure 
over momentum is d^p/(27r)^. For finite T, in the imaginary time formalism, the integral 
over four momentum goes over to an integral over three momentum and a discrete sum 
in the manner given in [3]. We will assume the reader is familiar with finite temperature 
field theory. 

We now return to problem of understanding what it means to use the one-loop cor- 
rected mass to calculate one-loop corrections. If we iterate the one-loop result once we 
obtain 



V{cf>)^V{cf>)-^z y^IVln 



Ik 



- Ml 



(7) 



If we keep only the contributions from (7) when the didj both act on the same mass- 
squared matrix we obtain, by carrying out the differentiation and expanding the log. 



^ n=l 



1 

p n 



A 

k 6 



(8) 



Here the term inside the square brackets is understood as a matrix with two indices; the 
trace is over powers of this matrix. It may be verified that the RHS is the tree-level plus 
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contributions from the daisy graphs of (1), with the correct combinatorics except for the 
two loop graph (n=l) which is over counted. The correct factor for the two-loop graph is 
half the one appearing from this procedure [14] . In the large N limit if we keep only the 
leading order in N and account for the two loop difference, this will reproduce the result 
of Dolan and Jackiw [3] for the finite temperature daisy sum correction to the effective 
mass at the origin: 

(9) 

Following [3] , one keeps only the most infrared divergent contributions. This amounts to 
letting the d/d(j? act only on [p^ — ^(0^ — f^)]""" for n > 2. For the two loop case it may 
act on either [p^ — ^(0^ — v"^)]^^ or [k"^ — ^(0^ — v"^)]^^ which again explains why this 
case is different. As shown in [3], however, the daisy graphs do not exhaust all dominant 
N contributions. For this one needs to sum the superdaisy graphs. Therefore, even 
for large N, this approximation scheme is not complete. Furthermore, in (7) there are 
additional corrections when the di and dj both act on different mass-squared matrices. 
These corrections cannot be written as a sum of IPI diagrams. Therefore a sensible 
approximation scheme can only be achieved by ignoring these corrections (which should 
in any case be subleading in N). In the last section we will compare our expansion 
result with the leading 0{N) correction from using an effective temperature dependent 
mass. 



2. (j)^ theory to subleading order. 

The systematic 1/N expansion allows us to calculate (3) as a perturbation in 
near = [15,16]. Root [16] has evaluated the leading and next-to-leading corrections 
to the zero temperature scalar potential in 4,3,2 and 1 dimensions. The procedure at 
finite temperature is very similar. 

At very high temperature, our results for four dimensions should be essentially those 
of a three dimensional euclidean field theory with a dimensionful 0*^ coupling. In three 
dimensions and zero temperature, the leading 0{N) potential has long been known [15]. 
It has exactly the same form as the sum of finite temperature superdaisy graphs that were 
computed by Dolan and Jackiw [3] for a four dimensional 0^ theory. An important point 
in our approach is that the effect of introducing an auxiliary field a is to shift the mass 
term and as a result there are no infrared divergences in this formalism. 

To proceed, we first set A = f/N. The 1/A^ expansion assumes / is fixed as 
increases, not A. Then by introducing a dimension two auxiliary field a we rewrite the 
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action (2): 




The auxiliary field has eliminated the 0^ term; the original form of (3) is easily recovered 
by use of the equation of motion for a. 

To calculate the effective potential V{(f)) one proceeds as follows. First, using the 
background field method one computes the effective potential as a function of backgrounds 
of and a. Then, the background of a is eliminated by its equation of motion. 

The systematic expansion is performed by expanding the action Sf^, a] about 
backgrounds (p and x thus: 

^i^^^i + ^i^ + (11) 

The factors of have been inserted with hindsight. In particular, the rescaling of the 
background is necessary in order to ensure that the effective potential is renormalizable 
to each order in the 1/N expansion. This will be seen later in the explicit calculations. 

One now writes a formula similar to (3), with an integral over and a: 

r[0,x] = S[VN<P,x]-tnin J [D^][Da]exp [j, j - -0>^a + . 

(12) 

The factor of % is useful in order to compare results in the expansion to those from 
the usual perturbation theory in h\ for practical purposes we will take h = 1. All terms 
linear in the quantum fields have been discarded, and we have defined 

A.^— + (13) 

The integral over is gaussian and can (formally) be performed exactly to give the 
leading 0{N) contribution to the effective potential. The integral over over also pro- 
duces a terms which may be expanded in a power series. It is known that the o"^ terms 
give rise to the next-to-leading corrections in the expansion, which is all we are 
interested in. Furthermore, to calculate the effective potential we can take and x 
space-time constants, but not a. The integral over produces a kinetic term for a which 
must be carefully determined in order to evaluate the next-to-leading corrections. 

The 00" term in (|r^) can be rewritten as a term purely in a by shifting the variable of 
integration 0. This can easily be acheived to all orders in the expansion, but since 
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we are only interested in the leading and next-to-leading terms it suffices to write the 
argument of the exponential in (|12]) as 



2^N 



(14) 



Using 0*(x) = —i^ sJ(x) ^^P(i / 'Jk4>'') evaluated at Ji = 0, and then rescaling the 
source and field (f) hj we obtain for the integral over in (|T^) , 



ih In / [D(f)] exp 



-z/iTrlnAr.i 
2 



i^ln 



2 ^-^ 2v/iV 



(7 



evaluated at J = 0, where 



1 



2VA^ 



,(15) 



(16) 



Of course, eq. (|T5]) does nothing more than define a perturbative expansion with a 
scalar propogator A*-'. The first term in (^) is a one-loop result. It has been evaluated 
many times before and is given by [3] 



o 



9 -h / Trln 

2 J x,p 

= - [Vo + Vt. 

J X 



(17) 



Vq is the zero temperature result which is divergent and must be regulated. With a sharp 
momentum cutoff A, 



Vn 



N 
32^ 



Vt is the finite, temperature dependent, result: 



Vt = N 



X 



-r Infy/^^l + ^—y^ 



;i8) 



(19) 



We have dropped O{0^) terms and higher, and c ~ 5.41. This result is valid at high 
enough temperatures. In the electroweak model, it is known [8] that the high temperature 
expansion is well justified at temperatures relevant for the study of the phase transition. 

The 0(o"2) correction is a 1-loop contribution as well. Extracting the IPI part from 
(p^) we obtain to leading order 



-a 



-X {P + idf - X 



a. 



(20) 
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In deriving this result we took the functional derivatives in (^), fourier transformed 
the resulting (5-functions and dropped all total divergences. We further used prop- 
erties of the translation operators: exp{—ip ■ x)H{—idp) exp{ip ■ x) = H{x — idp) = 
exp{—idx ■ dp)H{x) exp{idx ■ dp). Dropping total divergences assumes that the integral is 
well regulated (in three dimensions it is finite). Eq. (|^) may be evaluated using Feynman 
parameters; the result is 



o 



h I dz f — — -a , =a 



a, 



cr 



Jx 2^7^(3 _ z{z - 1)(92 



X 1^71(3 



sm 



1 + 4x/92 



(21) 



which is 1//? times the corresponding zero temperature result in three dimensions [16]. 
The p integral was evaluated using /^[p^ — x\~'^ = — c^^/c^x^ /pln[p2 — x] and only the 
leading temperature dependent and field dependent correction was kept. The value of the 
arcsine is uniquely determined since for = the integral over z is unambiguous. 

The final result for (|l^) to the required order is thus 



r[0,x] = S[^/N<P,x] + ^^^TrlnAr.i + ri[0,x], 
where the next -to-leading contribution from the gaussian integral over a is 



(22) 



ri[0,x] = -i^Trln 



1 + 



247r/? 



sm 



'l + Ax/d\ 



(23) 



We now have all the formalism behind us. Before presenting explicit calculations it is 
useful to rederive the leading order "mass-gap" equation of [3] to see how the infrared 
problem is automatically avoided. 



To leading order in one drops the a terms and eq. (|T2|) reduces to 



A^,^ 3A^ , . ,9 9^ 1 , . 1 



N 



- ^x{4> 



v^) + Uh [ Trln[p2-;^] 



(24) 



For convenience, we have rescaled the vev, v"^ Nv"^ with respect to the last section. 
The equation of motion for x, 







(25) 
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gives the following equation 

X = k<P'-v')-th^ I (26) 
6 6 jp — X 

But X is also just —2dT/d(f?, which is the radiatively corrected "effective mass" to leading 
order and is the reason why ( p6|) is called the "mass-gap" equation. It is exactly the result 
derived in [3] by summing superdaisy graphs, each of which is infrared divergent. The 
finite temperature, T = result in four dimensions is [3,15,16] 

where the ". . ." refer to 

(c - i - ln[A2/?2i) + (28) 

^2 ^ ^ ^^967r2^ 967r2 ^ ' 

This result is the sum of the T dependent finite part and the T = divergent part. The 
divergences may be absorbed by renormalizing v"^ and / and by introducing an arbitrary 
mass scale M [15]. Then for / ~ 0(1), and for all reasonable values of M(3 we can forget 
the extra terms at high enough temperatures. 

To compare, the zero temperature result for the action (2) in three dimensions is [15] 

^ = L[ct)''-v^)-h^-^. (29) 
^ 6 ^ 6 47r ^ ' 

Eq. ([27| ) determines T2, the temperature when the effective mass vanishes at the 
origin, to leading order in A^. Putting x = we get {h = 1) 

4 = I2v\ (30) 
P2 

Furthermore, at the origin, ^Jx has a remarkably simple solution for T just above T2 [3] 

X^^fi-iV (31) 



3 \(3 P2 

In general, we find for y/x the leading result 

/ 



487r/3 



1 + (1202/32 _ 12^;2/52 + 1) 



(32) 



We have chosen the sign in the solution (|3^) of (^7]) so that ^/x is positive for T bigger 
than T2. Eq. (^) simplifies in various limits. At 02 = 0, for T just above T2 we obtain 
(ll). At 3847r202/32 > / and T just above T2, 

V^-^M- (33) 
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Finally, when T/T2 - 1 > //(647r2) (i.e. 327r2(l - 12v^(3^)/f > 1) then for all (f) we have 
the simplification 




T2)2 + 2(T-T2)T2) + 



(34) 



We now return to the computation of Fi in the high temperature limit. Some care must 
be taken to ensure that the final answer contains no temperature dependent divergences 
[17]. Also, to the order we are working in, it is sufficient to use the solution of (P7| ) for x 
[16]. We will drop all field independent divergences. 

We make the replacement —>■ —p^ and the trace becomes an integral over space-time 
and momentum. At finite T,ijp^ -P'^ I d^f /C^ny^En = -['^t^^P]'^ Iiridf)VWEn 
with —p"^ —>■ + AiT'^n'^T'^ in the integrand. The discrete sum over n runs over all integers. 
One then obtains, 



In 



h 



1 + 



hf/i2AiT/3) 



+00 

E 



sm 



-1 



+ 



(35) 



The integral is not infrared divergent, but is ultraviolet divergent. The ultraviolet 
divergences should be the same as those appearing if we did the momentum integral over 
four momenta. [This not quite the same as taking 13 ^ 00 here, because we use the 
high temperature result for the a propogator.] We must isolate these divergences before 
calculating the temperature dependent part in the high temperature limit. Furthermore, 
even for say the n = contribution we were not able to find an analytic expression for 
the integral. We can however evaluate it in some limits. [For the rest of this section, we 
delete overall integrals over space-time.] 



Consider the n 
\) becomes 



contribution. We then change variables: t = \j^xlv^- Equation 



ri(n = 0) 



(4x)^ /--rft 



1 + at tan ^ ( - 
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3x + 4 



(36) 



where we have defined a = hf /{2A'K(3^/Ax) and used the identity sin (1 + t 



-1/2 ^ 

tan~^ t~^. At 0^=0, a is small to temperatures just above T2 because from ( |3TD we have 
a < 1 as long as 1 — f3/f32 > //(327r^) at the origin. In what follows, we will assume small 
a and small enough (f)^/x- In (P^), as t ranges from to 00, t tan~^ ranges from to 
1. At t = 1, ttan^^t^^ = 7r/4. Therefore for small a, attan'^t"^ remains small for all 
values of t in the integral. Hence, for small a and small 0^ we can use the simplification 
ln[l + r] ^ r to simplify (^) or (p5|). 
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For the moment let us ignore the exphcit term in Fi. To calculate ( P^D for small 
a we note the following observation for the n 7^ contributions. If ^Jx < ttT then the 
arcsine possesses an expansion in inverse powers of the momentum. We have, 



lv/l+4x/(p^+4.^n^T^)J ^ 7r/2 _^ (-l)'(4x)'+^ 

The n = contribution unfortunately does not possess such a straightforward expansion 
for all p^. If > 4x then the expansion is as above, otherwise it is 



We then split up Fi into four parts: A, the contribution from n = for > 4%; 5, the 
contribution from n = for < 4%; C, the contribution from n 7^ for > 4%; the 
contribution from 7^ for < 4%. i? is found using the expansion (^Sf), whereas all the 
rest involve the expansion (0). We denote the contributions to A, C and D from the first 
term in (^) by a, c and d, respectively. We further define 6 to be a contribution similar to 
a but with the momentum integral up to < 4%. Then it can be seen a + b + c + d is field 
independent and can be dropped, i.e. we evaluate Fi = {A~ a) + B + {C — c) + {D — d) — b. 
For B we get 

(4Y)i ^ f"^ dt 1 /1\ fh^Y , , 



For b we get 
We also get 



487r/?2 7o (27r)3 2v^ 967r2/?2- 



(40) 



487r/?;e^^'+'" (2/ + 1) 

.41) 

247r/5 Il67r2 12/32/ 567r3/32' ^ ^ 



where is given in the appendix (for e = ^4%, x = 0). Furthermore it is straightforward 
to show that D — d\s subleading in temperature and may be dropped. 

The 02 term from the expansion of the log is more easily evaluated and contributes 

where 1\ is taken from the appendix (e = 0, x 7^ 0). If we keep all powers of 02^ 
ignore the arcsine term, we must also include 

fix* ^ (-/0V6x-)'(a - 5)!! 
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This series converges for all values of cfP' /x we will be interested in. 

Each of the 0{h?) divergent terms in the above equations is not separately renormal- 
izable, however their sum is. This is because for the next-to-leading terms we can use eq. 
(0) to express 0^ _ h^/Ai^f] = Qx/ f + - h{12/3^)-^ + 0{l/N). Since the tree-level 
potential contains a term proportional to v'^Xi 'we see that the above divergences in Fi 
can be absorbed by renormalizing v"^ (throughout we have neglected field independent 
constants). The dominant contributions to Fi are found by summing the above results; 
we obtain 




^ ^/^'^y/X ^ fix 



-^x 



247r/? 
^X^ 



-/0V6x)'(2/-5)!!//!, 



;i + /0V3x)^-l 



(44) 



Altogether, summing the various partial results we find for the dominant high temper- 
ature effective scalar potential from (|12D, V{(f)) = — F[0, x(0)], the following renormalized 
result up to a constant {h = 1): 



V(0) = N 



:X(0' 



2/ 



X 



+ {N + 2) 



X 



3 

X^ 



24/^2 i27r/3 



;i + /0V3x)^+iV-l • (45) 



where x is a solution of dV/dx = 0. We do not give details of the renormalization which 
for zero temperature and in 3 and 4 dimensions may be found in [15,16]. 

We will use these results in the last section of this paper, where we also discuss the 
range of 0^ and T for which these results are the most dominant. 



3. Discussion. 

Our main result is the leading order in N and next-to-leading order expression for 
the high temperature scalar potential, eq. (|45|). The tree level potential, with the nor- 
malizations of the last section is 

V{<P) = ^{<l>'-vr. (46) 

Eq. (45) includes loop corrections to this. We assumed P^/x ^ 1 aiid P^/x ^ //487r ^ 
//150 in obtaining the next-to-leading result. At T = T2 for the leading order result, this 
gives approximately 100^7/ > 0^ > fv^lOO. At 0^ = we require T - T2 > /T2/3OO. 

To study the nature of the phase transition we may look for zeros of dV/dip^. There 
is always one zero at the origin, 0* = for all i. Away from the origin we can look for 
zeros of dV/d<p'^ = {dV/dx){,dx/d<f) + dV/d<i? = dV/d<i?. Since to 0{N), dV/d<i? oc x, 
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the leading order result has one such zero away from the origin when T < T2 (which is 
the local minimum at temperatures below T2). Thus, the leading order result does not 
admit a first order phase transition, for which we need zeros away from the origin for 
temperatures T > T2. 

At next-to-leading order, the critical temperature T2 is modified from the leading 



order result. dV/dclP' = for V given by ( ^5|) still has a solution when x = 0. Writing 
out dV/ dx at 0^ = and setting x = immediately gives for the temperature T2 , 

1 + ^) = ^^v'. (47) 

T2 has been reduced from its leading order value. dV/d(j)'^ = has one further solution at 
the origin. dV/d(j)'^ = gives 

fv^=^(l + /0V3x)l (48) 

when 7^ 0. At (/)^ = this has the solution ^ = // {12N7rf3). Writing out dV/dx = 
and inserting the above expression for gives a second (slightly higher) temperature 
when the effective mass at the origin vanishes. To 0(1/A^) it is given by 

+ ^"j^/^""' ^ t2 = 12v^ (49) 

For iV = 4, / = 1 this is within 0.5 percent of T2. Between this temperature and T2 the 
effective mass at the origin is negative! However, since our approximations break down 
at P^/x less than about //487r such a phenomenon is not necessarily physical; here it is 
an artifact of our approximations. 

To ascertain the nature of the phase transition we have to study the shape of the 
potential near T = T2. If a second minimum degenerate in energy with the one at the 
origin does not appear by temperature Ti > (1 + 5/ x 10~^)T2 then we cannot reliably 
determine if the theory admits a first order phase transition. However, in this case, if the 
exact theory does admit a first order phase transition then it will be only a very weak 
one. 

The general solution of x(0) is complicated. As noted by Root [16], one can use the 
leading order result for x(0) in the next-to-leading potential ( |45| ) to find V{(f)) to 0(1) 
in the 1/A^ expansion. However, (^) can develop an imaginary part for temperatures 
below the leading order value for T2. Therefore, it is preferable to use the next-to-leading 
order solution for xi^P)- Our results in what follows are correct to 0{1/N), barring the 
other approximations we made. At T near T2 given by (|47|) we once again have (|33|) when 
0^ ^ f^/300. To see this, we write out the equation for x(0), 

6X ^ 2_ 2 1 + 21N_VX 
f 12/?2 47r/? 

1 + /0V3x)^ - i] + TTT^Wrli + /</>V3x)i (50) 



X 
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^2 and y/x = ^/0^/6 + S. Assuming that S is small one obtains 
{\/3 — l)/N] which can be ignored when 5 <^ J fcjy^/Q. We can improve 



In this put (3 - 

5 ^ -//487r/?[l 

on this estimate for ^Jx by noting that when ^/x ~ \J f4>^ 1^ the last two terms in (|50|) 
add to an amount which for large enough (say > 3) are much smaller than the third 
last term. Hence we have the approximate solution 



/ 



487r/? 



12t;2/52 + 1 + 2/A^) -1 



(51) 



, which at 
This result holds 



This should be valid at temperatures near and 0^ sufficiently far from the origin. 

To see if a first order phase transition is possible we must find the zeros of dV/dcfP' 
dV/d(lP'. There is one solution at ^Jx = and another given by solving (| 

T = T2 and X ~ /0V6 has a solution at ^ ^ t;/y^3//[2A^(A^ + 2)] 
if A^ is not too large (\/0^ is not too small). At A^ = 4 we have ~ \/Jv/A. We can 
do better numerically. We studied the case A^ = 4, / = 1. The approximate solution 
(^) suggests Ti = T2(l + 3.5 x 10"'^), with ~ 0.08f at the second minimum at 
T = Ti. The height of the barrier separating the two minima is at most of 0(2 x 10~'^v'^). 
At T = T2 the minimum is at ~ 0.15f , which is slightly smaller than our cruder 
estimate above. 

All the numbers above occur as our approximations break down so we cannot determine 
the exact nature of the phase transition. However, as argued above we may deduce from 
them that the exact model can admit only a very weak first order phase transition. To 
get a more accurate picture one should evaluate (35) in the limit P^/x ^ //487r. 

To compare these results with the effective mass insertion procedure discussed in 
the introduction, we note that the leading 0{N) one-loop corrections to ( ^Uj) give the 
(renormalized) potential [3] 



V = N 



4! 



24/52 



(/(0^-t;^)/6)^ 
127r/3 



(52) 



As indicated in the introduction, this is complex for 0^ < ^2 rjj^g leading 0{N) one-loop 
corrected mass is 

(53) 




6 V" ^ ' 12/?2 J ■ 
Using this to perform the one-loop corrections, one gets 

_ ^2 + l/12/?2)/6 (/(02 _ ^2 + l/12/?2)/6) 



V = N 



24/?2 



127r/3 



(54) 



The effective mass at the origin vanishes at T2 = 12v^, as for the leading 0{N) result 
from the 1/A^ expansion. However, there is a higher temperature at which this also occurs. 
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Namely, if we put T = T2+t then for small t the appropriate equation is easy to solve. The 
result is that the effective mass at the origin also vanishes at T ^ (1 + //144)T2. Between 
this temperature and T2 the effective mass at the origin is negative. Thus compared to 
our results from the expansion, we see that the effective mass insertion technique is 
not reliable at temperatures less than about 1 percent above T2 even in the large N limit. 

In conclusion, for large enough N the model can admit only a very weak first order 
phase transition. For = 4 it would be interesting to see how 0{1/N) contributions 

modify the results above. At = 4 the next-to-leading corrections are certainly smaller 
than the 0{N) corrections and we do not expect 0(1/A^) corrections to compete with 
the 0{N) or 0(1) contributions. However, it would be necessary to perform an explicit 
calculation to check this behef. Finally, in the context of the standard model, one needs 
to add gauge fields A to this model. In fact, the contributions from gauge loops are 
expected to dominate over those from scalar loops in the weak scalar self-coupling limit. 
The main problem in the case of a gauged model is the gauge fixing dependence of the 
computations. Otherwise, in principle, is should be possible to determine the dominant 
contributions to the finite temperature scalar potential. For example, for SU (2), the gauge 
kinetic term has two types of interaction and [^l^Af]^, where i = 1,2,3) as 

well as an A^ interaction. Loops from the second A'^ interaction and loops from the A^ 
interaction are subleading to those from the first A'^ interaction. If we drop these last two 
interactions then what remains can be treated very similarly to the pure (p* theory. One 
must then determine what gauge fixing invariant quantities can be extracted from such 
a computation. It is expected that, if carefully done, physically measurable quantities 
should be gauge fixing independent. Work along these lines is in progress. 
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Appendix. 

The finite temperature quantities we need to evaluate in this paper are all of the form 



y n=—Oi 



{2n: 



-[/ + 47rVr2 + x]-' 



(55) 



for / a positive integer. The discrete sum runs over all integer n. We need only the cases 
e = and x 7^ or e 7^ and x = 0. We assume < 1, but > //(487r) ^ //150. 

For / > 1 we can find the answer by differentiating the / = 1 result / — 1 times w.r.t. 
X- To evaluate the expressions we follow [3]. The sum over n is exphcitly performed to 
give 



(27r)3 



+ 



2Vp' + X VWTx (eMPVWTx) - 1) 



(56) 



The first term in this expression is the T — piece and contains all the ultraviolet 
divergences; the second term is the T dependent piece and is finite. Changing variables, 
the T dependent piece is 



1 



27r^P^Jpe + P'^x (exp( + P'^x) - l) ' 



(57) 



In the high temperature limit, and for e = 0, we have the leading result of [3], 



/i(e = 0) 



A2 



+ 



X 



167r2 87r2 



In 



A 



+ 



+ 



X 



12/32 AnP 167r2 



1 
2 



21n[/3Vx] 



(58) 



The log pieces combine to give a single term proportional to xln[/3A]. Such log-divergent 
terms in the final answer must be renormalizable. The coefficient of this term is of 
0(/3y^/(27r)) down from the preceeding 0(T) term, and in our approximation all 0{x) 
terms can be neglected in comparison to 0{^yx/P) terms. Furthermore, to properly 
renormalize such terms in the effective potential we should keep subleading corrections to 
the a propogator, eq. (21), which we discarded. Therefore it is consistent to drop all the 
terms proportional to x- This is what we do in the following. 

For e 7^ we subtract from the above the temperature dependent piece 



27r2/32 



z^dz 



z' + P'x 
e 



ex]){^yz'^ + P'^x) - 1 



z'^dz 



27r2/32 Jo ^2 + ^2^ 



yz^+p'x 



^ tan-i 



27r2/3 27r2/3 
and the zero temperature piece, 

z'^dz 



eVe^ + X 



47r2 Jo V^2"+X 



.\/X. 

eVe^ + X 
87r2 



87r2 



+ 



X 

87r2 



In 



V^^Tx + e 



\fX 



X 

87r2 



In 



\fX 



(59) 



(60) 
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to arrive at 



. . X A2 1 

/i(e^O) = ^ + 



167r2 12/32 27r2/? 27r2/3 



tan ^ 



TT 

2 



(61) 



In the limit x ^ this expression has an expansion in positive integer powers of x- 
Therefore, its /th derivative w.r.t. x evaluated at x = is well defined. By expanding 
the arctan in this limit we obtain for / a positive integer and x = 0; 



A2 1 

+ 



167r 



12/32 



^2,-21 



Si-i + 



27r2/3(2/ - 3) ■ 



(62) 



/; for e = and x 7^ is simply found by differentiating ( |6TD w.r.t. x the appropriate 
number of times. The answer is, for / > 1, 



hix + 0) 



X5-' (2/ -5)!! 
47r/? 2'-i(/ - 1)!' 



(63) 
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